Inelastic scattering rates in d-wave superconductors. 
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The inelastic scattering rates of quasiparticles in a two-dimensional d-wave superconductor, which arise from 
interactions with either acoustic phonons or other quasiparticles, are calculated within second order perturbation 
theory. We discover a strong enhancement of scattering with collinear momenta, brought about by the special 
kinematics of the two-dimensional fermions with Dirac-like spectrum near the nodes of the d-wave order parameter. 
In the case of a local instantaneous interparticle potential we find that either an RPA-type resummation of the 
perturbation series or an inclusion of non-linear corrections to the Dirac spectrum is called for in order to obtain 
a finite scattering rate in the limit ui/T — ► 0. In either way, we find drastic changes in the scattering rate, as 
compared to the naively expected cubic temperature dependence. 



Recent measurements of microwave [T] and 
thermal conductivity [2] in 90 K YBCO have 
probed the quasiparticle relaxation rates in a 
range of temperatures below T c for which inelas- 
tic scattering is believed to be important, and 
rates approximately proportional to T A and T 2 , 
respectively, were revealed. Interestingly enough, 
both behaviours deviate from the simple cubic 
rate which would follow from the naive Golden 
rule estimate for the el-el scattering rate in the 
presence of a linear density of quasiparticle states. 
The likelihood of this discrepancy being due to 
an intricate difference between the quasiparticle 
lifetime versus, generally different, charge and en- 
ergy relaxation rates motivates one to revisit the 
calculations based upon the golden rule. 

We first consider the scattering of quasipar- 
ticles off three-dimensional acoustic phonons at 
temperatures below the Debye temperature. Af- 
ter having integrated the second order perturba- 
tion theory result for the inelastic quasiparticle 
lifetime over the out-of-plane phonon momenta, 
one recovers the expected cubic rate 
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which is frequency independent in the limit u> <C 
T. However, the prefactor 1(1) obtained from 
the angular part of the planar momentum inte- 



gral becomes divergent when the speed of sound 
s approaches the value of the quasiparticle ve- 
locity V2 parallel to the Fermi surface (the latter 
is proportional to the maximum superconduct- 
ing gap and therefore much smaller than vf)- 
7(1) ~ In(«2/|«2 - s\). 

Although posing no real problem for the el-ph 
scattering rate, the above divergence reflects a 
potential danger in the case of el-el scattering of 
quasiparticles with equal velocities. 

To this end, we consider a local instantaneous 
interaction between either charge or spin densi- 
ties. In the leading approximation, one can ne- 
glect scattering processes between any pair of the 
neighbouring nodes because of the strong eccen- 
tricity of the equal energy contours in momentum 
space (vp/v2 ~ 14 ± 3) [3]. For the remaining 
intra-node and opposite-node processes scaling of 
the local momenta near the nodes by the corre- 
sponding velocities gives rise to the symmetrical 
Dirac-like quasiparticle spectrum E(k.) = |k|. 

Further analysis reveals, that amongst the 
terms corresponding to different combinations of 
coherence factors [4] and nodes, involved in the 
scattering processes, the most important ones are 
those which correspond to scattering off ther- 
mally excited quasiparticles with nearly parallel 
momenta. 
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It proves convenient to write the scattering rate 
in terms of the charge or spin susceptibility 
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valid for <C T. Eq. ([2]) is then integrated with 
the energy conserving delta function and the ap- 
propriate combination of distribution functions, 
and for the external momentum taken right on 
the node this yields the integral 
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which tends to diverge in the limit lo/T — > 0, un- 
like the omitted regular terms proportional to cj 3 
or ujT 2 . A mere substitution of u> by T would 
yield a cubic temperature dependence but con- 
ceal the divergence stemming from the momenta 
q w 0. This kinematical singularity, which re- 
flects enhancement of scattering between parti- 
cles with collinear momenta, has also been noted 
in the early studies of semi- metals [5] . 

In lattice models, the above divergence is read- 
ily cut off by non-linear corrections to the bare 
quasiparticle spectrum. In the case of parame- 
ters relevant for the problem of high temperature 
superconductors (t 1 ~ — 0.5i) we obtain 
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Nonetheless, in the case of a strictly linear spec- 
trum no physically meaningful result can be ob- 
tained without a summation of the entire pertur- 
bation series. As an attempt to carry out this 
procedure one can substitute the bare interaction 
A by an effective interaction 
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and in general there might be other regimes of 
higher u> leading to different temperature depen- 
dences. Observe that neither this rate nor ([4]) 
reverts to a cubic temperature dependence when 
substituting T for uj. 

An important exclusion from the above generic 
situation is provided by the case of el-el inter- 
actions mediated by antiferromagnetic spin fluc- 
tuations, whose ordering vector is commensurate 
with the distance between the opposite pairs of 
nodes. In this case one arrives at a non-singular 



second order result r 



1/2 rp5/2 



Although it is widely believed that this partic- 
ular channel is relevant for the cuprates, it re- 
mains to be seen whether the lowest order esti- 
mate captures all the relevant physics. A formally 
related example of the system of interacting Dirac 
fcrmions describing quantum Hall plateau transi- 
tions 6J indicates that it may indeed be the case. 
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containing the singular susceptibility ^j. In the 
limit uj/T <C (XT/tA ) 2 , this procedure yields 



